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ABSTRACT: In this paper we generalize the quantum mechanical model of 
three coupled oscillators because of its usage in the crystal optical activity. Using 
this model we can include the influence of all essential couplings between single 
oscillators which represent the molecules or atoms of optically active crystals 
belonging for example to the space groups of symmetry D| or D%. The single 
oscillators are damped and therefore we can include both parts of the optical 
activity - optical rotatory dispersion and circular dichroism - into computations. 
We present more universal formulas for description of the above mentioned parts. 



1 Introduction 

In this paper we want to generalize the results of the previous paper |I[ in 
which we solved the optical rotatory dispersion (ORD) of crystals belonging to 
the space groups of symmetry D% or D%. The theory was based on the specific 
quantum mechanical model of three coupled linear harmonic oscillators represent- 
ing the elementary cell of these crystals as the smallest structure system. But 
the model used in [ED] was solved for the undamped single oscillators forming the 
compound oscillator and therefore this model described ORD only as one part 
of the optical activity (OA) and the results were valid only for the wavelength 
range outside the absorption band. Now we can use the model of three damped 
coupled oscillators and in this way we also describe the second part of the OA 
- circular dichroism (CD) - and remove the limitation of the wavelength range. 
But as will be shown, we need the results of the paper fTJ for the solution this 
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generalized model. 

ORD is the dependence of the rotation of the linear polarized light per unit 
length on the frequency or the wavelength. CD is the ellipticity of the outgoing 
wave from the optically active medium. It is based on the fact that the right 
and the left circularly polarized waves, in which the linear polarized wave is 
decomposed, and which propagate in optically active medium with different phase 
velocities, are absorbed differently. OA is fully described by the complex rotatory 
power p that is given by the relation 

p = p + ia = 7^(ni- n r ) = ^ [(nj - n r ) n r )] , (1) 

where p is the ORD and a is the CD, n\ and n r are the complex refractive indices 
of medium for the left and the right circularly polarized waves, ni, n r and Ki, n r 
are the real and the imaginary parts of complex refractive indices. 

It is well known that the molecules or atoms forming the crystals of the space 
groups of symmetry D\ and lie on the helix with screw axis of symmetry 
and their optical activity is caused by interactions among molecules or atoms in 
the most important cases. This is true (with the exception of camphor) for all 
typical representative crystals belonging to these space groups. In the paper [|IJ 
we gave reasons for using the three coupled oscilators model. Here we mention 
again that this model fully represents the elementary cell of the crystals belonging 
to the space groups of symmetry D| and and it also enables to include the 
couplings between even and odd oscillators on the helix. In we have proved 
that these couplings also have the substantial influence on the OA in the case 
when the diameter of the helix is comparable with the height of the elementary 
cell or smaller than its height. But it is known that it is valid for all important 
crystals with the screw axis of symmetry. 

Our model of three coupled oscillators is the same as in |]J. The oscillators 
lie on the helix, which axis (the crystal axis c) is parallel with the z coordinate 
axis. The second oscillator lies in the plane z — 0, its coordinates are x, y, 
and its direction of vibrations is described by direction cosines cos a, cos/5, cos 7. 
The first oscillator lies on the helix in the plane z = —d and its direction of 
vibrations is turned by the angle —6 around the z axis with respect to the second 
one. Similarly the third oscillator lies on the helix in the plane z = d and its 
direction of vibrations is turned by the angle +8 around z axis with respect to 
the second one. The angle 9 is of course 120 deg for the crystals belonging to 
the space groups of symmetry and D®. We suppose that all oscillators are 
identical and that we can neglect the influence of other helices which is given by 
the character of couplings between adjacent oscillators on one helix and between 
oscillators belonging to the different helices. Further we study the OA only in 
the direction of the crystal axis where OA is not screened by birefringence. 

Further it should be noted that the problems of crystalline OA (or ORD or 
CD) were also solved by using other methods. Here we can mention for example 
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the exciton theory [[3], f|, ||, |6|], Lagrangian formalism 0, Schwinger scattering 
H etc. But we think that the application of coupled oscillators models in the 
crystal optical activity is still important because the results are convenient for 
good approximation of the OA experimental data. This is also the most important 
aim of this paper. We want to get a new extension of the formulas for description 
ORD and CD experimental data. 

2 Crystal OA based on the three coupled oscil- 
lators model 

We can set the Hamiltonian for one compound oscillator which contains three 
coupled oscillators in the field of the left and the right circularly polarized waves 

-, rT fi 2 » 3 2 * 2T 
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+ \Mi (nr 2 + r 2 r 3 ) + /^r^] * + ]T F^p^ ^, (2) 

where ri,r 2 ,r 3 are the displacements of the oscillators from their equilibria, 
A i i( r i r 2 + r 2 r 3) + A i 2 r i r 3 is the potential energy of mutual interactions of the 
oscillators and F^ ,r = eE^ r ; £ = 1,2,3 are the electric forces projections of the 
left and the right circularly polarized wave into the vibration directions of the 
oscillators, e is the charge of electron, E l ^ r is the intensity of the electric field. The 
upper indices /, r hold for the left and the right circularly polarized waves, are 
the momentum operators of all single oscillators. The small positive parameter 
7o gives the possibility of the adiabatic interaction at the time t = —oo. The 
damping following from the limited lifetime of oscillators in their excited states 
is formally introduced by the parameter 70 . 

The electric vectors E l,r of the left and the right circularly polarized waves 
propagating along z axis have the components 

E l x = E e- i{ujt - kiz) , E l y = E e- i ( b *- kt *-%) , 

E r x = E e~ i{u)t - krZ \ E r y = E Q e- i ^ krZ+ ^) (3) 

and the vector E l ' r projections of these waves into the directions of vibrations of 
all oscillators are 



E l f = E e~ iuJt 



I 7T 



a cos 9 + (3 sin 9) e~ t<pl - r + (-a sin 9 + (3 cos 9) e l 
E 1 / = (a + Pe"*) , (4) 

a cos 9 - sin 9) e^" + (a sin 9 + (3 cos 9) e*^"^) 



E 1 / = E e~ iuJt 
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where the values 4>i,r = ki >r = n^ r ujd/c are complex phase shifts for the left (index 
I) and for the right (index r) circularly polarized waves. In all terms in eq. @ 
and further the + sign holds for the left and the — sign for the right circularly 
polarizad waves. 

Due to the mutual interaction the natural frequency ujq splits into three ad- 
jacent frequencies. We can solve our problem of three coupled oscillators using 
the normal coordinates qi,q2, Q3 which we can easy derive for the system of three 
classical coupled oscillators. If we have derived the relations for normal coordi- 
nates we can express the ri,r2,rs as combinations of these normal coordinates 
that are 

1 1 1 

ri = ; = gi + - 1 =q2 + -j==q 3 , 
l2 + A\ V2 J2 + Al 



r 2 = 1 qi + 1 93, (5) 

1 1 1 

r 3 = ; qi 7=F<?2 + , ?3 



2 + A\ V2 J2 + A 2 3 



where 



* = 2JS ' = ^ ' (6) 



After substituting (|5|) into (0) we get 

H * = -^? 1 ^ + -^? 1 9 i*-2( A4li4a ^ + A4a ^ + /4l>il ^)* 



2 + Aj y J V2 



+ p= (Fj' r + A 3 F^ + F^ 

; 2 + Al V 



(7) 



It may be simply verified that the expressions — y^^i^F^+y/^F^+F^), ^(F 1 /,r - 
F3' r ) and -, 1 = (F-[' r — \/2F2' r + F^) in the last equation are the electric forces 

in three normal coordinates. We denote them as F^ r , F^ r and F^ r . The equation 
(0) can be separated now by means 

Q2, qs, t) = t)tf 2 (g 2 , testes, *), (8) 
H l * = H\: + H\l + H l £. (9) 
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Substituting (§) and (|9|) into (0) we get three equations that are 

^q n ^r,{q v ,t) - ~ — 7^ + ^-q v ^ v {q v ,t) 

+— (F^e^^t), 17 = 1,2,3 (10) 

where 





= <4 + 


V2 + 




3/if 




2m 






= ul- 


^2 

m ' 






Wo 


= (4 + 


A*2 - 


V^2 + 


S/if 



2m 



We see that the natural frequencies of single oscillators ujq are split into three 
frequencies of the normal modes of vibrations. 

If we use the eq. (f|) we can derive for the electric forces in the normal 
coordinates the expression 

F% = e (ajf) Eoe-i^n) (12) 
and the coefficients a^ r are given by the relations 

2 q,2 I q2 

U l q r \ = — ZJL. L4 2 + 4A 1 cosfl + 4cos 2 fl±40 iir .sinfl(-,4 1 -2 cos 6 
" ' 2 + Ai 



(a l q l) 2 = 2(a 2 + /3 2 ) (sin#±20, ir cosfl)sinfl, (13) 

2 Q, 2 _|_ /J 2 

(V' 3 r ) = j^- [v43 + 4A 3 cos6» + 4cos 2 6'±40 Zir .sin6'(-A3-2cos6 

2 + v4 3 

cr^ - only have the meaning of the phase shifts. We can see that the Hamiltonian 
has a nonperturbed part 

^-~2m"^f + ~ 9 " (14) 

which depends only on the index rj of normal mode and does not depend on 
the polarization of the light wave. On the other hand the perturbed part of the 
Hamiltonian 

< = £j «) E p qv e-<^)^ (15) 
depends also on the polarization of the wave. 
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The mean value of the electric dipol moment induced by the left and the right 
circularly polarized wave that we hold as a small perturbation we can solve by 
means of the Kubo theorem ||. We get the relation for the mean value of the 
dipol moment 



K) *> = (e «) *>„ + £ «) E„ {{« n , K )l e-(»-S) + ,.., (is) 



but the first term on the right side of this equation is the constant dipole moment 
of the system and it has no meaning for us. The expression (^(q v ,p qtl )) is the 

Fourier transform of the retarded Green function ^qn,P qri )) of operators q v and 
p qv , that is 



n 1 e 



iut—^fot 



(17) 



For the retarded Green function \ iq n ->Pq n )) m the ground states of quantum 
system |^o) we hold the relation 

Jv>P*)) t = (vo | | Vo) • (18) 

In the eq. flT8|) is the unit step function, qJt) is the operator of q v in the 
interaction representation 



m° n tin* -m° tin 



q„(t) = W'V 



(19) 



and for the application of the operator in the exponent on the wave function 
holds 



D ±iH° q t/h 



e ±iuj Vn ti 



Vn)- 



(20) 



In following rearrangement we use that for matrix elements of the production 
of any operators F q and K q holds in the normal modes 



Vn 



F K 



V, 



Vk [Vk 



K. 



In 



(21) 



and further we use that the matrix elements of operator p q are 



Vn) = ™, m „ (r] m \q v \Vn) ; 



(22) 



UJri — UJn — U)ri 
'hnn 'Im 'In 



Now we can solve that 



,Vo 



Vo) = im 



(23) 
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and after substituting the eq. ( p3|) into ([18]) and solving the integral on the right 
side of eq. ([[7]) we have 

\\ _ 2ium ^ u} Vko \(Vk\q v \Vo)\ 2 
We can solve the relation for the induced dipole moments dq r by means of 



D and flT|) 

d *,r = £ V^L/ ^l^'^''/ /' £uC -<^+^)+td* (25) 



2e 2 ( 




2 


^l^o)| 2 









k 

and we can introduce the oscillator strengths of the normal modes of vibrations 
into ( P5| ) by the relation 

_ grmj^ |(^|g^ho)| 2 

Now the induced dipole moments can be expressed as 

^W^A^, (27 ) 

We can get the mean polarizability x l q r °f the crystal volume unit taking the 
part of the induced dipole moment (the part at electric field) and multiplying it 
by effective number N' of compound oscillators in a volume unit. If we suppose 
that the last single oscillator in any compound oscillator is the first oscillator in 
the next compound oscillator on any helix then the effective number N' = N/2 
where N is the number of single oscillators in a volume unit @ . If we introduce 
X 1 ^ into Drude - Sellmaier dispersion relation for refractive indices of the crystals 

3 

^2 _ i _ A^^TsM 
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l,r x ^" Z_^i Ag^ 

v =l 



we get 



n 



l,r 



_ 2nNe 2 * , , r , 2 

If we use the eqs. (29) and (p73f ) we can solve the difference of quadrats of 
complex refractive indices for the left and the right circularly polarized waves 



_ 2 _ 2 8vriVe 2 



n t - n r 



*E 

k 



(a 2 + (3 2 ) sin0(02 + r ) 

n 2 cos e+A 3 r 
COS ft n i a'2 JnM.. 

(30) 



m 

2cos0+Ai £ f Q 2 cos g+ A3 x 

2+Af Jll k0 Jq2 k0 COS V_ 2+A^ ■/93 fc0 
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and if we further use the relation (j^) for the complex rotatory power in which we 
use (pi + (f) r = ujd{fii — n r )/c and nf — n 2 — (fii + n r )(fii — n r ) we get 



AnNde 2 



-uo 



(a 2 + (3 2 ) sin#-]T 



k L 



2 cos 9+ Ax r 
2+Aj J1i k0 



2i^qU) 



2 cos 9+M f 
2+Ai JQ3 k0 



ul„ - oo 2 - 2i-f u 



OJ 2 — H^qLO 



(31) 



Now we may find the real and imaginary part of p(uo) that are from ([I]) ORD 
and CD 



Kep(u) = p(uo) = Auj 2 ■ 



2 cos 0- 



2+A 



-Af n (oo 2 -oo 2 ) 



\9(l 



/<Bu„ COSy [ U 2- - 



k I 

2 



2,. ,2 



2 cos 
24 



e+A 3 f ( 

^Al~J<K k0 V 



(32) 



lmp(cu) = (j(uj) = 2A"fou £ 



E 

k 



2cos6>+Ai 

J i 



2+Af 



in 



UO 2 ) + A 1 2 oo 2 



2 cos fl+ A 3 J? 
2+A^ Jl3 k0 



Ko-^ 2 ) + 4 ^ 2 H„-^ 2 ) +47o 2 ^ 2 



(33) 



where A 



AnNde 2 
mc 2 



(a 2 + [3 2 )sm6. 



We know that in the case \uq — uo 2 \ 2 > ^qOO 2 (it means in the frequency 
region far from dichroic frequencies) the ORD result (p2|) must pass into ORD 
result in the paper [Q where we used the model of three undamped coupled 
oscillators. The ORD result in the |L| is expressed in the form of three terms with 
defined frequency dependencies. The first term has the frequency dependence 
of the type uo 2 /(oJq — oo 2 ) and it is known as the Drude term. The second term 



contains dependence uo 2 / (uOq- 



-oo 



2\2 



(the Chandrasekhar term) and the third one is 



proportional to u 2 / (oo 2 — u 2 ) 3 . These terms frequency dependences are multiplied 
by the following coefficients 



K 



K, 



(i) 

fcO 

(2) 
kO 



2cos9 + A 1 2cos9 + A 3 

hi kn + Jl2 kn COSt> o , A2 ■/«W 



Pi 

m 



2 + Af 
2 cos 9 + Ax p2 
2 + A 2 A3jq ^o + ^Jvh 



2 cos 9 + A 3 



2 + A 2 lJq:i ™ 



2 + A 2 
cos 9 



(34) 
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K, 



(3) 
fcQ 



\ii 2 cos 9 + A\ 



m 2 + Aj 

H 2 2 cos 9 -\- A% 



A ifgi,„ -2/g 2 , cos 9 



A*i 2 + A 2 

We can rewrite our results fl32|) and ([33|) in similar way and we get more conve- 
nient results for the approximation of the ORD and CD experimental data 



^ 2 -E 

k 



'(1) 
L fcO 



Kir! [00 



2 

kO 



K 



(^o-^ 2 ) Z + 47o 2 ^ 2 



+ 



(2) 
AO 



K^kO ~ U 



2\2 









1 




3 

> 



(35) 



<r(o;) 



2A 7o cu 3 • £ 

k 



+ 



/v 



(1) 
fcO 



u; 2 ) 2 + 4 7 ^u; 2 



+ - 



2^? (u& - c 2 ) 



+ 4 7o 2 a; 2 



K {3) 

^feO 


3 (^ 2 - ^) 2 - 4 7o 2 ^ 




Ko-^ 2 ) 2 + 47o 2 ^ 


3 



(36) 



where we have neglected the expressions that are multiplied by small values [if, 
fif, a4 and also the expressions that are multiplied by p\ and p\ in the 
numerators of the first terms in compound brackets. 



3 Discussion 

We can see that our results (|3l|), (|32]) and ( |3~3|) do not contain the rotational 
strengths. But the nonzero values of rotational strengths are the condition for 
the existence of the optically active quantum transition and the results of OA 
are often expressed by using these strengths (see ref. [TC]). We are going to show 
that our results can also be expressed in this way. 

The rotational strengths are defined as the imaginary part of the scalar prod- 
uct of the electric and the magnetic transition dipole moments that means 

R k0 = Im ((ofd\k) ■ (k\tfi\0}) (37) 

where |0) is the ground state, \k) any excited state and d, in are the operators of 
the electric and the magnetic dipole moments. We will compute the rotational 
strengths in the normal coordinates in our case of three coupled oscillators 



R iv k0 = Im ( (Vo\dg v \Vk) ■ (Vk\m qri \r}o) ) , 77 = 1,2,3. (3 
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We have solved the normal components of the rotational strengths for the 
model of three coupled oscillators in ffl. We have derived the expressions 



R 



in 



R 



R 



he 2 d (a 2 + f3 2 ) sin 9 (2 cos 9 + A x ) 
2mc (2 + ~Aj) 
he 2 d (a 2 + j3 2 ) f q2kQ sin 9 cos 9 
2mc 

Tie 2 d {a 2 + (3 2 ) f q3m sin 9 (2 cos 9 + A 3 ) 
2mc (2 + If) 



(39) 



and we can see that the complex rotatory power ([31]) can be expressed as 

Rn 



p{uS) 



he 



r)=l k 



'/AO 



— uj 2 — 2i^ u 



(40) 



In the same way we can also express the ORD and CD results 

ul -uj 2 )r„ 



p{uj) 



8ttNuj 2 Jt, _ 

he ^ ^ ( • 



'/AO 



qr lkO 



167riVur7o 
he 



ri—l k 
3„, 3 



'/AO 



— UJ' 



+ 47o^ 



2,. ,2 



EE 



ft 



»,-J /•■ (uj 2 m -uj 2 ) +A 1 2 uj 2 



(41) 



(42) 



and we can accept these formulas as the crystal analogs of the formulas which are 
well known from the optical activity of molecules [11| or were recently derived 
for delocalized molecular aggregates [O. 

The ORD and CD are often expressed in the dependence on wavelength in 
practice. Because of this reason the ORD and CD formulas that are rewritten 
in this dependence are also often used. We can present several ORD and CD 
formulas that are suitable for theoretical approximation of experimental data 
and that are based on the single terms of the formulas ([35]) and (^). For ORD 
we get the formulas 



P(A) 
P(A) 

P(A) 



E 



tfif') t \ 2 \ 2 



lk 



(A 2 - A 



kOJ 



~k (A 2 - A| ) 2 + T 2 A 2 



E 



(A 2 - A 2 ) 2 - r 2 A 



^2 \2 



(A 2 - A 2 ) 



+ r 2 A 2 



(43) 



E 



Km A 4 (A 2 - A| ) 


/A 2 - A 2 ) 2 - 3r 2 A 2 ; 




;(a 2 - a 2 ) 2 + w 


3 



where K\ p k \ K^ k , K^} are the constants defined only by the crystal parameters. 
The central dichroic wavelengths A^o correspond to the transition frequencies Uko 
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and the constant T contains the damping constant 70 of single oscillators. It is 
evident that forms of all these constants can be easy derived from fl35l) . The 
values of these constants are computed by the approximation of experimental 
data by using the least squares method. 

We can see that the first formula in ([|3]) passes into well known Drude for- 
mula for wavelengths very different from central dichroic wavelength. The second 
formula passes into Chandrasekhar formula [H| on the same conditions. It is 
also known that for example the ORD of a-quartz is very well approximated by 
the sum of the Drude formula and the Chandrasekhar formula with one central 
dichroic wavelength. The form of the first two relations in (|43|) was derived in 
14[] by means of the classical model of two coupled oscillators but the form of 



constants 

K (p) 

was different. 

We are able to determine similar conclusions also for the relations by means of 
which we can approximate the CD experimental data. The terms in the relation 
(|36|) can be rewritten into dependence on wavelengths by using new constants as 



a (A) 
a (A) 

a(X) 



k (A 2 - A 2 r + ha 2 



E 
E 



j^f) \ / \ 2 \ 2 



L 2k 



A 3 (A 2 - A 



kOJ 



[(A 2 - xio) 2 + r 2 A 2 



(44) 



E 

k 



^?A 5 


3(A 2 -A 2 ) 2 -T 2 A 2 ; 




/A 2 - A 2 ) 2 + r 2 A 2 ~ 


3 



The form of the first and the second relation was also derived in [13] with ex- 
ception of the form of all constant K^ 7 ' . We have made the first calculations which 
should tell us the meaning of these practical relations and above all the meaning 
of the third relations in the eqs. (fEJ) and (0). Now we can mention that the 
third relation in (|43D has a very little meaning in the approximation of the exper- 
imental ORD values of a-quartz. If we approximate the ORD values by means of 
all three relations (for one central dichroic wavelength) we can essentially reduce 
the sum of least squares and better approximate the ORD data mathematically 
but the value of computed central dichroic wavelength is not probable. On the 
other hand if we approximate the CD experimental data of sodium uranyl acetate 
crystal (SUA) JT^| by temperature 77 K we get much better results by using all 
three relations (|44D with comparison to using only the first two relations. The 
experimental data are very well described by theoretical formula and the values 
of all constants are physically real. But the practical meaning of relations ( |43|) 



and (54) can be verified in the future by the approximation of ORD and CD 
experimental data of more crystals. 
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